Nematic Valley Ordering in Quantum Hall Systems 
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The interplay between quantum Hall ordering and spontaneously broken "internal" symmetries in 
two-dimensional electron systems with spin or pseudospin degrees of freedom gives rise to a variety of 
interesting phenomena, including novel phases, phase transitions, and topological excitations. Here 
we develop a theory of broken-symmetry quantum Hall states, applicable to a class of multi- valley 
systems, where the symmetry at issue is a point group element that combines a spatial rotation 
with a permutation of valley indices. The anisotropy of the dispersion relation, generally present 
in such systems, favors states where all electrons reside in one of the valleys. In a clean system, 
the valley "pseudo-spin" ordering, or spatial nematic ordering, occurs via a finite temperature 
transition. In weakly disordered systems, domains of pseudo-spin polarization are formed, which 
prevents macroscopic valley and nematic ordering; however, the resulting state still asymptotically 
exhibits the QHE. We discuss the transport properties in the ordered and disordered regimes, and 
the relation of our results to recent experiments in AlAs. 



A remarkably diverse set of phases, exhibiting the 
Quantum Hall (QH) effect, are observed in sufficiently 
clean two-dimensional electron systems (2DES) sub- 
jected to a high magnetic field Of these, a partic- 
ularly interesting subset occurs in multi-component QH 
systems, where in addition to the orbital degree of free- 
dom within a Landau level (LL), electrons have low en- 
ergy "internal" degrees of freedom, such as spin or a 
"pseudo-spin" associated with a valley or layer index. 
QH states in such systems, in addition to the topologi- 
cal order, characteristic of all QH states, feature broken 
global spin/pseudospin symmetries - a phenomenon 
termed QH fcrromagnctism (QHFM). The entangling of 
the charge and spin/pseudospin degrees of freedom leads 
to novel phenomena in QHFM states, including charged 
skyrmions Q, finite temperature phase transitions |3j], 
and Josephson-like effects 0, Q . 

In the cases studied to date, the global symmetry is 
an internal symmetry that acts on spin/pseudospin. In 
this paper we study a situation where the global sym- 
metry acts simultaneously on the internal index and on 
the spatial degrees of freedom. This occurs naturally in 
a multi-valley system where different valleys are related 
by a discrete rotation, so that valley (pseudospin) and 
rotational symmetries are intertwined. An example of 
such a system which is central to this paper is the AlAs 
heterostructure 041], where two valleys with ellipsoidal 
Fermi surfaces are present, as illustrated in Fig. [If a). 

This linking of pseudospin and space in this system 
has two significant consequences at appropriate filling 
factors such as v = 1. First, in the clean limit the on- 
set of pseudo-spin ferromagnetism, which occurs via a 
finite temperature Ising transition, is necessarily accom- 
panied by the breaking of a rotational symmetry that cor- 
responds to nematic order, with attendant anisotropics 
in physical quantities. We shall call the resulting phase 



a quantum Hall Ising nematic (QHIN). (The Ising- type 
pseudospin fcrromagnctism is consistent with the general 
classification of anisotropies in QHFM Q). Second, any 
spatial disorder, e.g. random potentials or strains, neces- 
sarily induces a random field acting on the pseudospins 
which thus destroys the long ranged nematic order in the 
thermodynamic limit. Interestingly, though, the result- 
ing state still exhibits the QHE at weak disorder so we 
refer to it as the quantum Hall random field paramagnet 
(QHRFP). 

Although for concreteness we shall focus on the simple 
case of the v = 1 state in AlAs-heterostructures, our 
findings are readily extended to other values of v and a 
variety of multi- valley systems. 

Symmetries: The only exact symmetries of QH sys- 
tems are the discrete translational and point group sym- 
metries of the underlying crystalline heterostructures. 
However, in many circumstances, there are additional ap- 
proximate symmetries, some of which are continuous. To 
the extent that spin-orbit coupling can be ignored, there 
is an approximate U(l) spin-rotation symmetry about 
the direction of the magnetic field. Since the magnetic 

length, Ib = y|§ > and the Fermi wave-length, X F , are 
long compared to the lattice constant, the effective mass 
approximation is always quite accurate, so it is possible 
to treat the translation symmetry as continuous. If the 
electrons occupy only a valley or valleys centered on the T 
point in the Brillouin zone, the effective mass approxima- 
tion also elevates a C n point-group symmetry to a contin- 
uous U(l) rotational symmetry. Terms which break this 
symmetry explicitly down to the discrete subgroup come 
from corrections to the effective mass approximation, and 
so are smaller in proportion to (a/A^?) 2 , where a is the 
lattice constant of the semiconductor. All three of these 
approximate symmetries hold in GaAs heterostructures. 
However, once there are multiple valleys centered 
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on distinct symmetry-related Bloch wave-vectors, the 
effective mass tensor for each valley is, gencrically, 
anisotropic. Thus, already in the effective mass approx- 
imation, individual valleys do not exhibit full rotational 
invariance; there are only the original discrete set of ro- 
tations which are associated with a simultaneous inter- 
change of valleys. These discrete symmetries are unbro- 
ken for weak interactions in zero magnetic field. How- 
ever, we show in this paper that in the presence of a 
strong magnetic field they are spontaneously broken at 
certain filling factors. 

Specifically, in the two valley case considered explic- 
itly here (Fig. la) the Hamiltonian has an approximate 
Zi x U(l) invariance: The Z2 represents the operation 
of a 7r/2 rotation combined with valley interchange. The 
U(l) reflects an approximate conservation of the valley 
index, which is violated only by the exponentially small 
Coulomb matrix elements, Vi V , which involve the inter- 
valley scattering of a pair of electrons. The QHFM should 
thus exhibit a finite-temperature Z2 or Ising symmetry 
breaking phase transition, accompanied by a spontaneous 
breaking of the rotational symmetry from C4 to C2, i.e. 
to Ising-ncmatic ordering. It is important to note that 
although Vi V breaks the approximate U(l) symmetry, the 
Ising symmetry is exact. 

We should note that there is a well understood counter- 
example to our general argument concerning the lack of 
continuous symmetries in multi-valley systems which is 
realized at a (110) surface in Si. Here the 2DEG occupies 
two valleys centered at k = ±Q/2, Q being shorter than 
the smallest reciprocal lattice vector. In this case, the 
only rotational symmetry is a symmetry under rotation 
by 7r. Yet, to the extent the intervalley scattering, Vi V , 
can be neglected, this problem was shown by Rasolt et 
al to have an SU(2) pseudo-spin symmetry. This de- 
rives from the fact that, in this case, the effective mass 
tensors in the two valleys are identical. While in the case 
of interest to us there is only a discrete Z2 symmetry, 
due to the effective mass anisotropy in each valley, in 
the limit of small anisotropy there is a reference SU(2) 
symmetry which is only weakly broken. For clarity, and 
without loss of generality, we will in places consider this 
analytically tractable limit, although in reality, the mass 
anisotropy in AlAs is not small. 

Ising anisotropy: The single-particle Hamiltonian in 
each of the valleys, labeled by the index k = 1, 2, is given 
by H K = V ^- K ;^ /c) \ where K x = (K Q ,0) 
and K2 = (0, Kq) are the positions of the two val- 
leys in the Brillouin zone. We work in Landau gauge, 
A = (0, — Bx), in which eigenstates can be labeled by 
their momentum p y that translates into the guiding cen- 
ter position X = p y l 2 B ■ The lowest LL eigenfunctions in 
the two valleys are given by, 
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FIG. 1: (a) Our model band structure. Ellipses represent 
lines of constant energy in the fc-space There are two non- 
equivalent anisotropic valleys, 1 and 2. (b,c) Schematic rep- 
resentation of two types of order in the QHFM. The ellipses 
here represent LL orbitals in real space. 
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FIG. 2: Easy-axis anisotropy of the QHFM as a function of 
underlying mass anisotropy. 



where A 2 = {mi yX /m\ ty } = (m2 :V /m2 : x) is the mass 
anisotropy in terms of which u\ — l/«2 = A. 

Proceeding to the effects of the Coulomb interactions, 
we notice that the terms in the Hamiltonian that involve 
inter- valley-scattering processes require large momentum 
transfer, of order ir/a, and therefore they are small in 
proportion to a/ £3- In accord with that, we write the 
Hamiltonian as follows, 

h = H + H iv ,H = ^^nqK«(qK'«'(-q), (2) 



where S = L x L y is the system's area, p KK is the density 
component within valley k , I^(q) = 2 -j|- is the matrix el- 
ement of the Coulomb interaction, and Hi V denotes inter- 
valley scattering terms f30| , which we neglect for now. 

To account for the spatial structure of LL wavefunc- 
tions, we follow the standard procedure of projecting the 
density operators onto the lowest LL (see, e.g., Ref. 0]): 



(1) 



P KK (q) = F KK (q)p KK (q), F KK (q) = e 



(3) 
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where the magnetic translation operator is given by, 



p KK (<i) = Y,e iqa cl x+ c K , x _, X ± =X±^ 



<1>, 



In the limit of vanishing mass anisotropy, A — > 1, the 
Hamiltonian Hq is S'[/(2)-symmetric, so at filling factor 
v = 1 there is a family of degenerate fully pseudo-spin 
polarized ground states, favored by the exchange inter- 
actions. 

= Y[(ac[ x + P4 iX )\0), \a\ 2 + \(3 2 \ = 1. (4) 
x 

In this notation, the components of the nematic order 
parameter are given by n x = aft* + a* (3, n y = ia/3* — 
ia*(3 7 n z = \a\ 2 — |/3| 2 , where n 2 = 1. We can use the 
states ((4]) to obtain a variational estimate of the energy 
per electron of the system for different (uniform) values 
of the order parameter which should be reliable at least 
for A near 1. The result is 
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where 



Dx = (Ci + C 2 )/2, D 2 = (Ci - C 2 )/2, 



2 g(^ 
Ci j= 

7T \/A 



1/A 2 ) 



,C 2 = 



2A 



I + A 2 ' 



(5) 



(6) 



(7) 



K being the complete elliptic integral of the first kind. 
Clearly, when A ^ 1, the SU{2) symmetry is broken 
down to Z 2 x U(l) and the resulting QHIN indeed has 
an Ising (easy-axis) symmetry. The magnitude of the 
anisotropic part of the energy, D2A0, is pictured in Figj2] 
For the experimentally relevant case, A 2 w 5, k ps 10, 
the anisotropy reaches a relatively large value of 5 K at 
B = 10 T. Let us also note, for subsequent use, that the 
Ising symmetry can be explicitly broken in experiments 
by the convenient application of a uniaxial strain Q, 
which then acts as a valley Zecman field. 

Thermal properties: In order to understand the be- 
havior of the system more generally, and in particular to 
describe the properties of domain walls and excitations, 
we need to account for spatially varying order param- 
eter configurations. The classical energy functional for 
smooth configurations of the order parameter can be ob- 
tained approximately for |A — 1| <C 1 by the method of 
Ref. Q: 



E[n(v)j = f/ d 2 r(Wn) 2 - | J d 



2 2 



(8) 



where a 
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^-(A — 1) . The symmetric part of the 



0(A— 1). In writing Eq.©, we have neglected anisotropic 
stiffness terms of the form, ^ J d 2 rn z ((d x n) 2 — (<9j,n) 2 ), 
£f- Jd 2 r[3(V?i z ) 2 - (Vn) 2 ]. While these terms are also 
quadratic in the gradient expansion, in the limit |A— 1| <C 
1, the first term is at most cubic [9( in A — 1, such that 
PA ~ o((A — l) 2 ), while the second term is quadratic, 
PA' = °((A ~ and so they are much smaller than the 
gradient term we have kept. 

The nematic ordering temperature can readily be es- 
timated from Eq.®, which is precisely the continuum 
limit of the 2D Hciscnbcrg ferromagnet with weak Ising 
anisotropy. Consequently, T c vanishes for a = 0, but 
only logarithmically, due to the exponential growth of 
correlations in the Heisenberg model, 



ksTc ~ 47rp s log 1 [p s /a£ 2 B ] . 



(9) 



stiffness coefficient in Eq. © is given by, p s = - 



6V2tte£b 



Since, in reality, the anisotropy is not small, a more ro- 
bust estimate is just T c ~ p s . This puts it in the range of 
several Kelvin, well above typical temperatures at which 
quantum Hall experiments are carried out, which range 
from a few tens to a few hundred milli-Kelvins Q ■ 
Domain walls and quasiparticles: The topological 
defects of an Ising ferromagnet are domain walls, in this 
case domain walls across which the valley polarization 
changes sign. We obtain a domain wall solution by mini- 
mizing the classical energy in Eq. ([8]) to obtain the length 
scale Ln = which characterizes the domain-wall 

width, and the surface tension, J ~ ^/p s a, its creation 
energy per unit length. The domain wall solution ob- 
tained in this way spontaneously breaks the approximate 
?7(1) symmetry as the energy is independent of the choice 
of the axis of rotation of n in the plane perpendicular to 
n z . Naturally, since the domain wall is a one dimen- 
sional object, thermal or quantum fluctuations restore 
the symmetry, but at T = 0, and in the absence of ex- 
plicit symmetry breaking perturbations, what remains is 
a gaplcss "almost Goldstonc mode" and power-law corre- 
lations along the domain wall. A small gap in the spec- 
trum and an exponential fall-off of correlations beyond 
a distance £j„ are induced when the effects of the weak 
intervalley scattering terms, V% v , are included. 

In AlAs, the anisotropy is A 2 ps 5, so Lq is only 
30% greater than £g, which indicates that our treatment 
should be supplemented by microscopic calculations that 
can better handle a strong Ising anisotropy A varia- 
tional ansatz for a domain- wall can be constructed of the 
same form as in Eq. [4] by treating a and /3 as (complex) 
functions of X, with asymptotic forms (a, 0) — > (1, 0) as 
X —00 and (a, /3) — > (0, 1) as X — > 00. In the limit of 
large A, the optimal such state consists of a discontinuous 
jump between these two limiting values across the do- 
main wall, so that the domain wall width is simply equal 
to Ib- When inter- valley scattering is absent, such a 
wall, being a boundary between two different QH liquids 
(one with v\ = 0, v 2 = 1, the other with v\ = 1, v 2 = 0), 
supports two counter-propagating chiral gaplcss modes — 
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one with pseudo-spin "up" and the other with pseudo- 
spin "down." Coulomb interactions between the two 
modes turn this into a type of Luttinger liquid. This 
connects smoothly to the description obtained above in 
the limit of weak anisotropy, and indeed the Luttinger 
liquid action can be derived explicitly from a er-model 
description [13, II | by augmenting the classical energy in 



Eq. ([8]) with an appropriate quantum dynamics. 

The other excitations of interest are charged quasi- 
particles and it is well known that in the SU(2) limit 
at A = 1 they are pseudospin skyrmions of divergent 
size 0]. However, the smallncss of Lq at A 2 ~ 5 alluded 
to above implies that for the experimentally relevant case 
the quasiparticlcs will be highly, if not completely, valley 
polarized. 

Properties of the clean system: For T < T c , where 
the pseudospin component n z has a nonzero expectation 
value, C4 rotation symmetry is spontaneously broken to 
CV Thus, non-zero values of any non-trivial traceless 
symmetric tensor can also be used as an order parameter. 

Ideally, thermodynamic quantities, for instance of the 
difference in the valley occupancies, provide the concep- 
tually simplest measures of the broken symmetry. How- 
ever, such quantities are not easily measured in prac- 
tice. Following our remark above, we should just as well 
be able to use the experimentally accessible transport 
anisotropy ratio 



N 
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^0 



(10) 
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as a measure of nematic order. However, at T = 0, where 
<j aa = 0, N is ill-defined. This problem can be resolved 
by measuring a aa at finite temperature, T > 0, and then 
taking the limit T — > 0. (Alternatively, one could imagine 
working at finite frequency, and then taking the limit 
as the frequency tends to 0.) However, in practice, the 
conductivity is strongly affected by the presence of even 
weak disorder, so any practical discussion of the resistive 
anisotropy must be preceded by an analysis of the effects 
of disorder. 

Length scales from weak disorder: By analogy 
with the random-field Ising model [l2| , we know [lH that 
even an arbitrarily weak random valley- Zeeman field de- 
stroys the ordering of the QHIN, leading to formation of 
"Imry-Ma" domains of opposite valley-polarization. In 
AlAs such disorder can stem from random strains, which 
lead to position-dependent relative shifts of the energies 
of the two valleys. While the average strain (i.e. the 
average pseudo-magnetic field) can be externally con- 
trolled [4J, fluctuations of the strain are inevitable. Ran- 
dom fluctuations of the electric potential, V, also give 
rise to a random valley field. 

The coupling of random strain and potential disorder 
to the QHIN order parameter is 

Est = d 2 rh(r)n z (r), (11) 



where the random field h(r) = [h st (r) + h pot (r)] with 
h st (r) tx — q , u(r) being the displacement of 

point r of the crystaf, and 

{»> 



h P ot(r) = 



m yY\ 



2irh 2 



dV 

dx 



dV 
dy 



(12) 



On the basis of this analysis, we expect that the ran- 
dom valley-Zeeman field is smooth, with a typical cor- 
relation length £dis ^ £b- For weak disorder, the Imry- 
Ma domain size is set by the mean-squared strength of 
the random field (assumed to have zero mean and to be 
short-range correlated) W = J d 2 r(h(r)h(0)) and the do- 
main wall energy per unit length, J (defined above) as 
£jm °c exp[ A(J) 2 /W } where A is a number of order 
1. Because of the exponential dependence on disorder, it 
is possible for £jm to vary, depending on sample details, 
from microscopic to macroscopic length scales. 

Disorder also leads to scattering between the valleys al- 
though this is again suppressed due to the mismatch be- 
tween the reciprocal lattice vector and the length scale of 
the dominant potential fluctuations. There is thus a sec- 
ond emergent length scale, which is the length scale 
beyond which conservation of valley pseudo-spin density 
breaks down. However, this length scale approaches a 
finite value in the limit of vanishing disorder due to in- 
tcrvallcy Coulomb scattering, discussed above. Different 
regimes of physics are possible depending on the ratio of 
£im/£iv Finally, especially when the filling factor de- 
viates slightly from v = 1, there is a length scale, £qp, 
which characterizes the quasi-particlc localization length. 
Because the magnetic field quenches the quasiparticle ki- 
netic energy, even for extremely weak disorder, we expect 
£qp ~ is relatively short. 

Intrinsic resistive anisotropy: In a quantum Hall 
state at low temperatures, dissipative transport is usu- 
ally due to hopping of quasiparticlcs between localized 
states, accompanied by energy transfer to other degrees 
of freedom [l4[ . Typically, transport is of variable-rangc- 
hopping (VRH) type, such that the optimal hop is deter- 
mined by the competition between the energy offset of 
the two states and their overlap. We will now apply these 
ideas to our system when its transport primarily involves 
hopping of electrons between localized states within one 
of the valleys. This requires a) either that a uniaxial 
strain be applied to substantially eliminate domain walls 
and achieve valley polarization in the proximity of v = 1 
or that the sample be smaller than £jm and b) that £™ 
be large compared to £,qp- For each valley the localiza- 
tion length is anisotropic, owing to the mass anisotropy, 
which results in the anisotropy of the corresponding con- 
tribution to the VRH conductivity. 

The contribution to the resistive anisotropy from 
quasiparticles in valley 1, N\, can be computed as fol- 
lows: First, we transform the anisotropic VRH prob- 
lem into the isotropic one by rescaling coordinates, x = 
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x/y/X, y = VXy. In the new coordinates the effective 
mass tensor is isotropic, which, given the uncorrelatcd 
nature of the potential, implies that the VRH problem is 
isotropic [3l| . and therefore <r xx = a yy . Since the ratio 
of the conductivities in the original coordinates is given 

J (Tyy 1? Syy ' 

* - (is) 

which is negative for A > 1, as expected, i.e. it is 
more difficult for particles to move in the direction of 
larger mass. Clearly, the resistive anisotropy produced 
by quasiparticles in valley 2 is N% = — Ni. 

At v = 1 localized states in both valleys are present, 
and due to combined particle-hole/valley-reversal sym- 
metry of the state (in the absence of Landau-level mix- 
ing), the density of localized states should be same: The 
resistivity is thus expected to be isotropic! However, for 
v ^ 1, particle hole symmetry is broken. Consider the 
case in which n z = +1, which corresponds to filling val- 
ley 1 states. Then, at slightly filling factor, v = 1 — 8v 
with 1 ^> 5v > 0, the density of localized states for val- 
ley k = 1 exceeds that for valley k = 2. Due to the 
exponential sensitivity of the VRH conductivity to the 
density of states, this implies that the contribution of 
valley 1 to the total conductivity dominates, leading to 
an anisotropy of the total conductivity N ps Ni. Con- 
versely, for v = 1 + 8v, N ps N2 — —Nil It is worth not- 
ing that the scaling argument presented above for VRH 
regime is likely more general, and also applicable to the 
regime of thermally activated transport, which is relevant 
at intermediate temperatures. 

Domain walls and the QHRFPM: Wc now move 
away from the above limit to where domain walls are 
a significant contributor to the transport — to systems 
much bigger than £jm and at weak uniaxial strain. Now, 
dissipative transport is complicated by the existence of 
multiple emergent length scales. Transport within a ne- 
matic domain proceeds by variable range hopping and/or 
thermal activation of quasiparticles. For length scales 
larger than m , it is likely to be dominated by transport 
along domain walls, which will have insulating character 
or metallic character depending on whether viewed at 
distances large or small compared to £ iv . 

A key question is whether the QHE survives the for- 
mation of domains. This is trickiest when no net val- 
ley Zccman field is applied where in the thermodynamic 
limit the domain walls form a percolating network. In the 
limit f;u — > 00, the associated edge channels are conduct- 
ing, and the domain wall network can be expected to be 
well described by two copies of the Chalker-Coddington 
network model [15[ at criticality. This implies a critical 
metallic longitudinal DC conductivity of order e 2 /h and 
the absence of QHE. However, at length scales longer 
than t;i V (or temperatures less than Ti V ~ the 
domain wall states are localized, which implies a phase 



that exhibits the QHE without Ising/Nematic order — the 
QHRFPM. Needless to say, in the absence of substantial 
amount of short ranged disorder (which can produce a 
relatively small value of the topological (quantum 
Hall) order in the QHRFP is likely to be fragile . 

(In some ways similar results were obtained by Rap- 
sch et al. [lj|, who considered an SU(2)-symmetric dis- 
ordered QHFM, where magnetic order is destroyed by 
forming a spin glass without destroying the QHE. ) 

In the presence of a uniform valley Zeeman field h the 
existence of the QHE is much more robust. Even weak 
fields can restore a substantial degree of valley density po- 
larization as domains aligned with this field grow while 
those aligned opposite to it shrink. Consequently, the 
domain walls no longer percolate but rather are sepa- 
rated by a finite distance that grows with increasing h. 
While we have yet to construct a detailed theory of the 
transport in this regime, it is clear that the characteristic 
energy scale characterizing the dissipative transport will 
rise rapidly from Tj„ for h = 0, to the clean-limit gap 
Ao ~ p s for substantial values of h. It is also important 
to note that this equilibrium response will come embed- 
ded in a matrix of dynamical phenomena characteristic 
of the random field Ising model that can be translated 
straightforwardly to the case of the Ising-nematic, as has 
been discussed in another context in Ref . (l7| . In partic- 
ular, the macroscopic nematicity induced by the appli- 
cation of h will be metastable for long times, even upon 
setting (h) = — thus giving rise to hysteresis. 

Experiments: Turning briefly to experiments we note 
that an anomalously strong strain induced enhancement 
of the apparent activation energy at v = 1 has been 
observed[4j in AlAs, where it was tentatively attributed 
to the occurrence of valley skyrmions. As we noted ear- 
lier, in view of our estimate of a large Ising anisotropy 
skyrmions of the requisite size (about 15 flipped pseu- 
dospins) are implausible. We would like to suggest that 
it is more plausible that these remarkable observations 
are associated with the growth of QHIN domains. In sup- 
port of this idea, we have estimated the domain size from 
the long ranged part of the potential disorder alone and 
find that it should be order the distance to the dopant 
layer, £dis ~ 50 nm which is thus much smaller than the 
system size. However, we currently lack a plausible es- 
timate of £j„ which is sensitive to the short ranged part 
of the disorder and which is needed to round out this ex- 
planation. Direct measurements of resistive anisotropics, 
and of hyseretic effects characteristic of the random field 
Ising model [ijj could directly confirm this proposal. 
Related work: We note that there is a sizeable body of 
existing work on Ising QHFMs produced at level cross- 
ings of different orbital LLs, which is typically achieved 
by applying tilted magnetic fields. These systems exhibit 
enhanced dissipation at coincidence jl8l421j which is the 
analog of a dissipation peak at zero valley Zeeman field 
in our language. Qualitatively, our results are consistent 
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with this earlier work. Where we differ is in our con- 
tention that the domain walls do not, even at zero valley 
Zeeman field, produce dissipation at T = — in the previ- 
this was not explicitly addressed in part 



ous work [22J, [22 
as the focus was on accounting for the unexpected dissi- 
pation. The reader will also note that the QHIN stud- 
ied here differs from the "nematic quantum Hall metal" 
(NQHM) phase which has been observed 24j in ultra- 
clean GaAs-GaAlAs heterostructures for fields at which 
the n > I st Landau level is nearly half filled. Unlike our 
system, the NQHM is a metallic state which does not ex- 
hibit the QHE, but has a strongly anisotropic resistivity 
tensor. 

In closing: The distinctive feature of our system is the 
breaking of a global symmetry that combines spatial and 
internal degrees of freedom. This physics and its atten- 
dant consequences will generalize immediately to other 
ferromagnetic fillings in the present system and then to 
other experimentally established examples of multi-v alley 
systems such as monolayer and bilayer graphene 25|, |26| , 
where two valleys are present, and Si (111) [13] , where, 
depending on the parallel field, either 4 or 6 degener- 
ate valleys can be present. Potentially, our ideas could 
apply farther afield in the case of 3D Bi, where three elec- 
tron pockets related by 2tt/3, Att/3 rotations are present. 
Recently, high-field anomalies in transport and thermo- 



dynamic properties of Bi were found |28|, l29( , which may 
indicate spontaneous breaking of the Z3 valley symmetry 
driven by magnetic field, reminiscent of QHFM. 
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